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Why are X-ray sources in the M31 Bulge so close to Planetary 
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ABSTRACT 

We compare a deep (37 ks) Chandra ACIS-S image of the M31 bulge to deep 
[0 III] Local Group Survey data of the same region. Through precision image 
alignment using globular cluster X-ray sources, we are able to improve constraints 
on possible optical/X-ray associations suggested by previous surveys. Our im- 
age registration allows us to rule out several emission-line objects, previously 
suggested to be the optical counterparts of X-ray sources, as true counterparts. 
At the same time, we find six X-ray sources peculiarly close to strong [0 III] 



line emission, the chances of these six X-ray sources lying so close to cataloged 



> 

emission-line sources, classified as PNe by previous optical surveys. Our study 
■ shows that, while the X-rays are not coming from the same gas as the optical 

O" 

PNe is only ~1%, suggesting that there is some connection between these [0 III] 
^ ■ emitters (possibly PNe) and the X-ray sources. We discuss the possibility that 

^ these nebulae are misidentified supernova remnants, and we rule out the possi- 

bility that the X-ray sources are ejected X-ray binaries. There is a possibility 
that some cases involve a PN and an LMXB that occupy the same undetected 
star cluster. Beyond this unconfirmed possibility, and the statistically unlikely 



one that the associations are spatial coincidences, we are unable to explain these 
[O III] /X-ray associations. 

Subject headings: galaxies: individual (M31) — planetary nebulae: general - 
X-rays: general — supernova remnants — astrometry 



1. Introduction 

The high-spatial resolution of the Chandra X-ray Observatory is allowing optical coun- 
terparts to be found for a large number of individual extragalactic X-ray sources. Such a 
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capability provides the first opportunity to find new examples of exotic sources for which 
there are only a handful of Galactic examples. Because of its proximity, low extinction, and 
high mass, M31 is an excellent place to search for such objects, and associations between 
X-ray sources and optical line-emitting sources provide interesting places to start. 

Early optical emission-line surveys of M31 (e.g. Rubin et al. 1972) and X-ray surveys 
of M31 (e.g. van Speybroeck et al. 1979) produced a limited number of reliable counterpart 
candidates due to resolution and sensitivity issues. More recently, many digital surveys have 
discovered hundreds emission-line sources, including several hundred compact [0 III] A5007 
sources in the M31 bulge, believed to be planetary nebulae (PNe) (Ford & Jacoby 1978; 
Ciardullo et al. 1989). 

At the same time, X-ray spatial resolution has improved dramatically, with Chandra pro- 
viding source locations to a few tenths of an arcsec. Because the emission-line sources in the 
M31 bulge are well separated, even in ground-based images, comparisons of high- resolution 
X-ray positions with optical emission-line images of the M31 bulge reveal interesting coun- 
terparts. Typically these sources of strong X-ray emission and optical line emission are found 
to be supernova remnants (SNRs) (Kong et al. 2002a, 2003; Williams et al. 2004), like the 
X-ray bright SNRs in our Galaxy. 

Several recent surveys of M31 with Chandra (Kong et al. 2002b; Williams et al. 2003) 
have revealed possible X-ray counterparts to optical emission-line sources classified as PNe 
by their appearance as point sources in [O III] images. Such counterparts are difficult to 
understand because PNe are among the weakest X-ray sources in the Galaxy. None of the 
Uhuru or HEAO A-l X-ray sources was ever optically identified as a PN. In fact, only in 
recent years has it been possible to study Galactic PNe using ROSAT and Chandra. Guerrero 
et al. (2000) detected a total of 13 PNe with ROSAT; all are faint and have very soft spectra. 
Chandra observations have provided more detailed information: NGC 6543 and NGC 7293 
have peak temperatures (MEKAL model) of 0.5 keV and 1.0 keV, and X-ray luminosities of 
10 30 erg s _1 and 3 xlO 29 erg s _1 , respectively (Guerrero et al. 2001). The most luminous 
planetary we are aware of is NGC 7027 with L x ~ 1.3 x 10 32 erg s -1 ; its temperature is 
~ 0.3 keV (Kastner et al. 2001). 

Another possibility is that these PNe are not counterparts at all, but instead are sepa- 
rated by a few arcsec, close enough to be confused with counterparts at the accuracy limit 
of simple catalog cross-correlation. To check this possibility, it is necessary to have images 
of the PNe sources to register and compare to the X-ray data. Recently, the Local Group 
Survey (LGS) team (Massey et al. 2001) has released their deep, wide- field emission-line 
images of M31. These images contain X-ray emitting globular clusters (GCs) as well as 
previously-cataloged PNe. These GCs allow precise registration of the [O III] and Chandra 
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images, reducing the X-ray error regions to a few tenths of an arcsecond. 

In this paper, we compare a deep Chandra ACIS-S image with the LGS [0 III] image 
of the M31 bulge in order to tightly constrain the association between several of the X- 
ray sources of the M31 bulge and bright, compact [O III] emitting regions classified as PNe. 
Section 2 describes the data and analysis techniques used. Section 3 discusses the comparison 
of the X-ray positions and the optical emission-line sources. Section 4 discusses possible 
explanations for the associations, and section 5 provides a summary of our conclusions. 

2. Data 

2.1. Data Processing 

We obtained the [O III], [S II], Ha, V and i?-band images of the M31 bulge from field 5 
from the Local Group Survey (LGS 2 ). These images have already been properly flat-fielded 
and the geometric distortions removed so that the coordinates in the images are good to 
~0.25" on the FK5 system, and the images through different filters are registered with one 
another. We therefore were easily able to subtract the V-band continuum from the [O III] 
image in order to make the [O III] sources stand out. The central 40" are saturated in the 
continuum image, rendering that section of the data useless for our purposes. Fortunately 
only 2 of the PNe counterpart candidates lie in this region. The [S II] and Ha images were 
used only to test the flux ratios as a diagnostic for photo-ionization of the X-ray PNe. 

We performed a rough calibration of the LGS [O III] image by matching the [O III] 
fluxes of 10 planetary nebulae (PNe) with published [O III] fluxes Ciardullo et al. (1989). 
This calibration provided a conversion factor of 5.5 xlO" 16 erg cm -2 ct -1 . We also roughly 
calibrated the Ha and [S II] images by matching the fluxes of the SNR DDB 1-15 (D'Odorico 
et al. 1980) to the fluxes measured in the calibrated data set of Williams et al. (1995) 
(Ha = 7.3 x 10~ 14 erg cm -2 s" 1 ; [S II] = 5.5 xlO -14 erg cm -2 s" 1 ). This calibration 
yielded conversion factors of 1.0 xlO" 16 erg cm" 2 ct" 1 and 1.8 xlO" 16 erg cm -2 ct -1 in Ha 
and [S II] respectively. Using these [O III], Ha, and [S II] factors, we converted the LGS 
count rates to units of erg cm" 2 s" 1 . 

We also obtained a deep Chandra ACIS-S image centered on the M31 nucleus. This 
image had an exposure time of 37.7 ks. We created exposure maps for this image using the 
CIAO script mergeall, and we found and measured positions for the sources in the image 
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using the CIAO task wavdetect. We detected 153 sources, 137 of which were located in regions 
outside of the saturated portion of the LGS images. The source list reached a flux limit of 
~2.5 x 10~ 7 ct cm~ 2 s" 1 (0.3-10 keV), or ~8 xl0~ 16 erg cm -2 s -1 assuming an absorbed 
power-law spectrum with slope 1.7 and Nh = 10 21 cm 2 , or a (unabsorbed) luminosity limit 
of ~7 xlO 34 erg s -1 in M31, assuming a distance of 780 kpc (Williams 2003). 

The X-ray image contained all but one of the previously mentioned PNe counterpart 
candidates (Kong et al. 2002b; Williams et al. 2003), which are all labeled and marked with 
their X-ray position error circles in Figure 2. The one counterpart candidate not detected 
(rl-23, Williams et al. 2003) is an X-ray transient candidate (Kong et al. 2002b). This source 
is in the saturated central region of the LGS data; even if it had been detected in the ACIS-S 
data, its position relative to the PNe could not have been further constrained by this study. 

2.2. Image Alignment 

We aligned the coordinate system of the ACIS-S image with the LGS coordinate system 
by translating and adjusting the plate scale of the ACIS-S coordinate system so that 13 
globular cluster sources had the same coordinates as the centroids of the respective globular 
clusters in the LGS V^-band image. This transformation, performed using the IRAF 3 task 
ccmap, had root-mean-square residuals of 0.16" in RA and 0.15" in DEC. These errors 
were added in quadrature to the position errors of the sources, determined by wavdetect, to 
calculate the final position errors for X-ray sources on the LGS [O III] image. 

Figure 1 shows the error circles of four of the globular cluster X-ray sources used to align 
the X-ray and optical images plotted on the [O III] image. The quality of the alignment is 
within the position errors of the globular cluster sources. The globular cluster X-ray sources 
in the aligned images are within la of the centers of the optical globular clusters. On the 
other hand, the PN sources, shown in Figure 2, are clearly not within the calculated errors 
of the X-ray sources. These slight offsets between the X-ray sources and the [O III] nebulae 
show that the [O III] emission is not coming from the same source as the X-ray emission, 
but the proximity of the sources in so many cases is intriguing, as it suggests the sources are 
connected in some way. 



3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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2.3. X-ray Spectra 

To further investigate the nature of these peculiar X-ray sources, we extracted their 
X-ray spectra from the ACIS-S observation using the CIAO task psextract. We binned the 
spectra so that they contained £ 10 counts per bin, allowing standard x 2 statistics. Source 
rl-2 contained sufficient counts to provide ^ 20 counts per bin. 

We then fit the background-subtracted spectra from 0.35-7 keV using the CIAO 3.0/Sherpa 
package (Freeman et al. 2001). We tried fits using Raymond-Smith (RS), blackbody, and 
power-law models, including absorption and a correction of the instrumental response for 
the well-known contamination build-up on the ACIS detectors. The results of the fits are 
listed in Table 1 and discussed below in §3.2. 

3. Results 

3.1. X-ray source locations 

We visually searched the aligned [O III] image for nearby X-ray sources by placing 2.5" 
radius circles onto the [O III] image centered on the locations of all X-ray sources detected 
in the ACIS-S images. This search yielded 15 counterpart candidates, 6 of which have been 
cataloged as PNe. Figure 2 shows these 6 "PNe" sources within 2.5" of the X-ray source 
position. The la error circles for the X-ray positions are plotted on the images. 

Five of the other 9 nearby [O III] sources are cataloged SNRs (r2-56, r2-57, r3-63, r3- 
69, and r3-84 (Kong et al. 2002a, 2003; Williams et al. 2004). The [O III] sources near 
CXOM31 J004229. 1+412857 (Kaaret 2002) and r3-87 (CXOM31 J004226. 1+412552, Kong 
et al. 2002b) are uncatalogued. The [O III] source near CXOM31 J004220.5+412640 (Kaaret 
2002) is a known emission line source (WB92a 26, Walterbos & Braun 1992) and radio source 
([B90] 64, Braun 1990). The source near rl-16 is a cataloged emission line object of unknown 
nature (Wirth et al. 1985). 

While all of these [O III]/X-ray associations are of interest, the SNRs have been studied 
in detail in previous publications, and the other sources have yet to be classified. Herein we 
focus on the 6 X-ray sources within 2.5" of "PNe". Of these objects, r3-67 is also within 2.5" 
of a known radio source ([B90] 122, Braun 1990). The radio source is more closely aligned 
with the X-ray source (1.7") than the [O III] nebula (3.7"). 

In order to determine the number of random X-ray/optical associations that would be 
found by searching within 2.5" of detected X-ray sources, we shifted the positions of our 
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X-ray sources by 15" and repeated the visual search for counterparts. This experiment 
was performed 4 times in 4 directions of position shifts. The mean number of random 
association found was 3, with an average number of cataloged PNe associations of 2. We 
therefore expect that 2 of our associations with PNe are chance superpositions; however, 
this number still suggests that 4 of the X-ray sources are somehow associated with PNe. 
According to Poisson statistics, there is only a 1.2% chance that all 6 of these associations 
are chance superpositions. 4 

For a more sophisticated determination of the likelihood of 6 false PN/X-ray associa- 
tions, we performed monte carlo tests for chance associations as a function of the separation 
between the PN and X-ray source. Stepping from maximum separation angles of 1" to 3" 
in 0.1" increments, we searched for PN/X-ray associations by comparing the X-ray and PN 
catalogs after rotating them about the M31 nucleus with respect to one another by a ran- 
dom angle (20-340 degrees). Any associations found in the rotated catalogs could not be 
real; therefore these searches provided the number of expected false associations for each 
maximum separation angle tested. 

We performed the monte carlo test using 2000 random rotation angles for each maximum 
separation angle. This number of trials provided a robust mean number of false associations, 
shown in Figure 3. In addition we monitored the percentage of tests that produced 4 or 
more false associations as well as the percentage that produced 6 or more false associations. 
These results are also shown in Figure 3. 

Our monte carlo tests confirmed what we found with our searches by eye, finding a 
mean of 2.07 false associations within 2.5" separation and 0.73 false associations within 1.5" 
separation. The monte carlo results also support our use of Poisson statistics for estimating 
the likelihood of false associations, as 1.3% of tests find 6 false associations within 2.5" and 
0.4% find 4 false associations within 1.5". Furthermore, as shown in Figure 3 we find 6 
associations at a range of maximum separation angles (2.3"-2.7"). Over this range, the 
probability of finding 6 false associations changes rather steeply (see Figure 3), covering a 
range of 0.4%-3.25%. This probability range represents our most conservative measurement 
of the probability that all of these associations are false. 

Finally, our precise image alignment provides position errors for the X-ray source lo- 



4 A 2.5" search radius was chosen as the root-sum-square of the typical position errors in the Chandra 
and Ford catalogs. It should be noted that if we run the same test with a 1.5" search radius, we find 4 PNe 
associations and an average of 1 random association. These numbers leave a 1.5% chance that all 4 of the 
associations within 1.5" are chance superpositions, suggesting that the statistics are not strongly dependent 
on the exact search radius chosen. 
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cations on the LGS [O III] image of only a few tenths of an arcsec. These small position 
error circles are shown in Figure 2. All of our X-ray sources are slightly, but significantly 
offset from the PNe positions. This offset is intriguing. It suggests that, while the X-ray 
emission is not likely to be coming from the PNe itself, the X-ray sources are likely somehow 
associated with the PNe. 



3.2. X-ray properties 

We analyzed the spectra of all six X-ray sources that lie within 2.5" of PNe. The spectral 
parameters for the absorbed power-law fits, including model type, model parameters with la 
errors, x 2 l v i an d absorption-corrected 0.3-7 keV luminosities, are given in Table 1. In one 
case, r3-7, the best fit absorption fell below the known foreground absorption. We therefore 
accepted the best fit model with the absorption fixed at the foreground value of 6 x 10 20 cm~ 2 
(e.g. Kaaret 2002; Trudolyubov et al. 2002). 

In one case, an absorbed RS model provided a fit as good as the fit from an absorbed 
power-law. Source r3-21 is as well fit by an absorbed RS model with Nh = 9±3 xlO 21 cm -2 
and kT = 4±2 keV. Even in this case, the absorbed power-law model is as likely as the RS 
model to be the correct interpretation of the data. Therefore, Table 1 shows the power-law 
fit, as it does for the other 5 sources. 

In the other 5 cases, an absorbed power-law model provided the only fits with x 2 l v < 
1.5. In three of these cases, the absorbed power-law fit was improved by the addition of a 
blackbody component to the spectrum. We applied an F-test to the x 2 values of the spectral 
fits with the addition of a blackbody component (e.g. Bevington & Robinson 2003). The 
F x values for rl-26, rl-2, and r3-67 were 23.1, 15.6, and 21.1, respectively (the others had 
F x < 1). These values leave <1% probability that a blackbody component should not be 
included in the fits for rl-26 and rl-2, and <5% probability that the component should not 
be included in the fit to r3-67. The blackbody components comprise 68%, 24%, and 99% of 
the modeled unabsorbed luminosities (0.3-7 keV) of rl-26, rl-2, and r3-67, respectively. 

The sources all show somewhat similar spectral properties. For example, all have power- 
law components with slopes consistent with 1.8-2.0, and all but r3-21 are consistent with 
Nh < 2 x 10 21 cm~ 2 . Even the absorption column of r3-21 is within 1.8a of this value. 

The X-ray variability of these sources has been investigated in Kong et al. (2002b). 
Object rl-24, rl-2, and r3-7 were found to have varying intensity in X-rays. The other 3 
sources have constant X-ray flux. These properties allow that these sources could be XRBs. 
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3.3. Optical properties 

We estimated the narrow-band fluxes of these sources in the LGS images. These fluxes 
were measured in 2.5" diameter apertures centered on the position of the [O III] nebula as 
determined by the IRAF task imcentroid. The counts in each nebula in each band were 
measured using the IRAF task phot and then converted to flux units with the calibration 
discussed in §2. The resulting fluxes are provided in Table 2. 

The [S II] /Ha ratio of line-emitting regions is a well-known diagnostic of shock heating, 
where ratios >0.4 are typical of shock heated regions, such as SNRs, and ratios <0.4 are 
typical of photo-ionized regions, such as H II regions or PNe (Levenson et al. 1995). Ford 322, 
Ford 21, Ford 13, and Ford 201 have [S II]/Ha ratios of 0.9, 1.4, 1.1, and 0.7 respectively, 
suggesting that these nebulae are shock-heated, and not photo-ionized. We are not able to 
determine the [S II] /Ha ratio for Ford 494 due to lack of a [S II] or Ha detection. Only for 
Ford 209 is the [S II] /Ha ratio consistent with a PNe hypothesis. 

The [O III] /Ha ratios (see Table 2) were corrected for absorption because of the large 
difference in wavelength. We estimated the extinction using the relationship of Predehl & 
Schmitt (1995). We assumed a standard reddening law, and we assumed Ay ~ A\p m\ and 
Aji ~ Atf a . We compared our [O III] fluxes and [O III]/Ha ratios to the objects in common 
with Ciardullo et al. (2002). Our [O III] flux estimates for Ford 13 and Ford 21 are in 
agreement with their measurements, as is our estimate of the [O III]/Ha ratio for Ford 13. 

On the other hand, our [O III] /Ha ratio estimate for Ford 21 is larger than theirs by 
a factor of 3. This suggests that our Ha flux estimate could be significantly low, causing 
the very high [S II] /Ha ratio (1.4). This discrepancy underscores the difficulty of measuring 
faint source fluxes in the high background of the M31 bulge. Ford 21 is the faintest PN we 
detect in Ha (1.1 x 10~ 15 erg s _1 ) in an area of very high background (background/source 
= 112). For comparison, Ford 209, the faintest PN in our sample in Ha (0.9 x 10~ 15 erg 
s _1 ), has a background/source flux ratio of 13. The high background of Ford 21 hinders 
the measurement of precise fluxes. We note that, although our flux ratio estimates are not 
precise enough to yield reliable source classifications, even if the Ha flux is truly a factor of 
3 stronger than our estimate, the [S II] /Ha ratio is ~0.5, still high for a PN. 

[O III] /Ha ratios are not reliable for indicating shock- heating or photo-ionization, as 
photo-ionized regions cover a wide range of [O III] emission strengths (Ho et al. 1993; Veilleux 
& Osterbrock 1987; McCall et al. 1985); however, it is interesting to note that four of the 
[O III] /Ha ratios (Ford 322, Ford 13, Ford 201 and Ford 209) are <5. These values are 
not out of the ordinary for interstellar shocks containing a variety of shock velocities ^ 100 
km s -1 (e.g. Vancura et al. 1992; Fesen et al. 1997; Mavromatakis et al. 2000). The other two 
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nebulae have [O III]/Ha ratios >10, suggesting photo-ionization, very high shock velocities, 
or O-rich emission regions. If these sources are SNRs, they could be young with high shock 
velocities: a characteristic consistent with their small sizes. 

While full optical spectra will be required to determine the true nature of these sources, 
the suggestion of shock heating from the narrow-band flux ratios provides further evidence 
that 4 or 5 of these 6 objects are not PNe. 

4. Discussion 

Because there is very little likelihood that these 6 X-ray sources lie so close to classified 
PNe by chance, we assume for the moment that at least some of these sources are somehow 
associated with the nearby [O III] nebulae. It should be noted that, even if the positions 
were coincident, there is no current model for, or example of, PNe that produces strong 
X-ray emission. There are also no models for X-ray binaries that produce strong [O III] 
emission. The unique Galactic X-ray source GX1+4 emits in [O III]; however, even this 
extreme example has [O III] /Ha <C 1 (Chakrabarty & Roche 1997). 

If these are true X-ray/ [O III] associations, what are they? The spectral indices and 
absorption indicated for the X-ray sources are typical of X-ray binaries. However, there are 
no PN formation models that would predict the presence of a nearby XRB. If these [O III] 
nebulae are SNRs, the XRBs may have been formed in the same supernova event, but we 
argue below that even this scenario is unlikely. 

4.1. Are the [O III] Nebulae SNRs? 

Before we address the nature of the X-ray sources, we attempt to understand the nature 
of the nearby [O III] nebulae. Their [S II] /Ha ratios and proximity to strong X-ray sources 
make these nebulae very strange indeed, as they are unresolved, strong [O III] sources in 
M31 that are likely shock- heated. In addition, any hypothetical X-ray sources coincident 
with the [O III] nebulae are not detected by Chandra, requiring that they are fainter than 
~10 35 erg s _1 in X-rays (see §2.1). 

There is a possibility, as several of the [S II] /Ha flux ratios suggest, that the [O III] 
nebulae are SNRs. If these [O III] nebulae are SNRs, they must be young. To be unresolved 
in the 1" seeing of the LGS images, these SNRs would need to be £4 pc across, smaller 
than the sizes of SN 1006 (~7 pc, Willingale et al. 1996; ~19 pc, Winkler et al. 2003) and 
3C 58, thought to be the remnant of the supernova of 1181 C.E. (~5 pc, van den Bergh 
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1990). These comparisons limit the ages of the SNRs to < 1000 yr. 

The low X-ray luminosities ( £ 10 35 erg s _1 ) do not rule out the possibility that these 
nebulae are SNRs, despite their implied ages. While the recent X-ray detections of several 
SNRs in M31 (Kong et al. 2003; Williams et al. 2004) show that some SNRs are located in 
regions of the ISM with densities >0.1 cm~ 3 , Magnier et al. (1997) found that most SNRs 
in M31 are weak in X-rays, suggesting that the density of the ISM is lower in the vicinity of 
many M31 SNRs. 

Furthermore, young SNRs of both Type la and Type II can be weak in X-rays. For 
example, SN 1006 and 3C 58 are thought to have been Type la (Winkler et al. 2003) and 
Type II (Panagia & Weiler 1980) supernovae, respectively. SN 1006 has an X-ray luminosity 
(1-10 keV) of ~5 xlO 34 erg s _1 (Winkler & Laird 1976; Winkler et al. 2003), and 3C 58 
has an X-ray luminosity (0.1-4 keV) of only ~10 34 erg s _1 (Helfand et al. 1995). Therefore 
neither of these young SNRs would have been detected by the Chandra surveys. 

Finally, if the [O III] nebulae are SNRs, some of them must have old progenitors in 
order to explain their location in the M31 bulge. Three of these nebulae (Ford 13, 21, 
322) lie within a 2' x 2' box centered on the nucleus, in the old stellar population of the 
M31 bulge. Although the M31 bulge is known to contain some molecular clouds (Melchior 
et al. 2000) and SNRs (Sjouwerman k Dickel 2001; Kong et al. 2003; Williams et al. 2004), 
allowing the possibility of some recent star formation, it is not known to harbor a young 
stellar population (Deharveng et al. 1982; Stephens et al. 2003). If these [O III] nebulae are 
SNRs, their progenitor stars were likely old. 

Because the M31 bulge does not contain a substantial number of young, high-mass stars, 
if these "PNe" are in fact SNRs, we must consider the possibility that some are SNRs from 
Type la supernovae. Such supernovae likely have old, white dwarf progenitors (e.g. Hoeflich 
et al. 1998). Type la supernovae do not arise from the core-collapse of a massive star; they 
do not produce stellar-mass, compact objects (i.e. XRBs). Even if these inner three [O III] 
nebulae are SNRs, the fact that the progenitor supernovae did not likely produce XRBs 
provides the first hint that these [O III] /X-ray associations cannot be SNR/XRB pairs. 

4.2. Are the X-ray Sources Ejected XRBs? 

Returning to the six strong Chandra X-ray sources displaced from the [O III] nebulae, 
we address the question of how these sources could be related to the neighboring nebulae. 
The X-ray sources certainly cannot be related to the nebulae if the nebulae are PNe, but 
if the nebulae are SNRs, the X-ray sources might be ejected X-ray binaries. To check this 



- 11 - 



hypothesis, we assume that the supernovae were able to create XRBs, temporarily ignoring 
the possibility that, if the nebulae are SNRs, some of them had Type la progenitors. 

The idea that these slightly misaligned sources are ejected binaries is not entirely radical. 
Such an association has been suggested in the Galaxy for the X-ray binary Cir X-l. While 
the initial suggestion of an association with the nearby SNR G321.9-0.3 has been ruled out 
(Mignani et al. 2002), there is a new suggestion that the radio synchrotron nebula which 
surrounds Cir X-l is actually its birth SNR (Clarkson et al. 2004). 

If these [O III] nebulae are SNRs, the comparisons of their sizes to SN 1006 and 3C 58 
(see §4.1) limit the ages of the SNRs to £ 1000 yr. If these SNRs produced the nearby 
XRBs, and the XRBs are 2.8-8.7 pc away from the SNRs' centers, these XRBs must be 
moving at least 2700-8500 km s -1 . Such velocities would be higher than those observed 
for Galactic neutron stars, which typically range from 100-1000 km s _1 (Cordes & Chernoff 
1998). 

Single neutron stars with very high velocities ( ^ 3000 km s _1 ) may conceivably exist. 
Arzoumanian et al. (2002) suggest that 15% of neutron stars have velocities >1000 km s _1 
and 10% of neutron stars younger than 20 kyr will lie outside of their birth SNR. In addition, 
selection effects in Galactic pulsar surveys could lead to an underestimate of the number 
of high- velocity neutron stars (Cordes & Chernoff 1998). Although such claims favor the 
existence of high-velocity neutron stars, there is no evidence for neutron star velocities >3000 
km s _1 . 

Moreover, it is difficult to imagine how a binary system could be born with a velocity of 
>3000 km s _1 . Any pre-supernova binary would have an orbital velocity that is only a small 
fraction of 3000 km s" 1 , and the asymmetric kick received by the compact object during the 
supernova event would be so fierce that it would unbind the binary. We therefore conclude 
that, whether the [O III] nebulae are SNRs or not, it is unlikely that even one of these six 
X-ray sources could be a high-velocity X-ray binary ejected from a SNR. 

4.3. Do the X-ray/PN Pairs Reside in Star Clusters? 

There is an additional non-random way in which a PN and an X-ray source could be 
so close to one another. The PN and XRB could be in the same star cluster, putting 
them spatially near one another while being unrelated physically. In such a scenario, the 
progenitors of the XRB and PN must have been formed at the same time. The neutron star 
in the XRB would have then formed from its high-mass progenitor before the PN formed 
from its lower mass progenitor. Applying this possibility to the M31 associations requires 
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that the star clusters themselves have gone undetected. 

There are Galactic examples of LMXB/PN associations. We cross-correlated the Galac- 
tic PN catalog of Kohoutek (2001) with the Galactic LMXB catalog of Ritter & Kolb (2003). 
This cross-correlation revealed only two associations with separations less than 250" (at typ- 
ical Galactic distances 250" corresponds to ~2.5" at the distance of M31). Both of these 
associations involve a PN and an LMXB that lie in the same globular cluster. The PN 
PK 065-27 1 and the LMXB 4U 2129+12 are 30" apart in NGC 7078 (M 15). The PN 
JaFu 2 and the LMXB 4U 1746-37 are 30" apart in NGC 6441. There is a second PN inside 
the tidal radius of NGC 6441 (7'; Bahcall & Hausman 1976); PN H 1-36 has a separation 
from 4U 1746-37 of 310". Monte carlo tests, in which we rotated the Galactic catalogs with 
respect to one another by random angles (0.1-5 degrees), show that in 1000 iterations there is 
<0.1% probability of 2 chance associations with separation less <100" and 5.5% probability 
of 2 chance associations with separation <250". 

Apparently, PNe and LMXBs that are close to one another in the Galaxy occupy the 
same globular cluster. However, globular clusters like NGC 7078 and NGC 6441 are easily 
detected in ground based images of M31. For example, NGC 6441 would have V 17.5 if 
located in M31, and would be easily detected like other globular clusters within V of the 
M31 nucleus with V £ 18 (e.g. ACH 6; Auriere et al. 1992). No such clusters are detected 
at the locations of any of the [O III] /X-ray associations. 

Clearly if these M31 LMXB/PN associations are located in clusters of stars, they are 
not globular clusters. In recent years, the distinction between globular clusters and open 
clusters has been blurred. Studies of extragalactic star clusters are revealing that they may 
form according to a power-law mass function (e.g. the Antennae; Zhang & Fall 1999). Each 
cluster then stays bound for a length of time determined by its mass, concentration, and 
orbital parameters (Fall & Zhang 2001). The fact that the star clusters in the Galaxy happen 
to fall neatly into two categories (open and globular) is due to the cluster formation history 
of the Galaxy. Other galaxies have different cluster formation histories and are likely to 
contain some star clusters that do not fit into these two categories. Examples of different 
cluster types in local galaxies include the "blue globular" clusters in the Large Magellanic 
Cloud (Hodge & Schommer 1984) and the massive and compact young star clusters in M31 
(Williams & Hodge 2001). 

The star clusters hosting the LMXB/PN associations in M31 would need to be the 
diffuse, faint relics of once-larger clusters. However, the likelihood that such relic clusters 
would be massive enough to contain both an LMXB and a PN simultaneously, massive 
enough not to dissociate completely in ~10 9 years, and faint enough to be undetected, 
seems low. 
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5. Conclusions 

Several recent Chandra surveys of M31 have identified X-ray sources in the M31 bulge 
spatially coincident with previously cataloged PNe. Using the [O III] data from the Local 
Group Survey, we have directly compared the X-ray source positions to the locations of the 
cataloged PNe, finding that the X-ray sources are not precisely coincident with the PNe. 
The estimated [S II] /Ha ratios of 4 of these "PNe" are consistent with shock- heating, i.e. 
SNRs. 

The fact that the X-ray emission is displaced from the [O III] emission is consistent with 
the classification of these nebulae as PNe, which are notoriously feeble X-ray sources. At the 
same time, the spatial proximity is unlikely to be coincidental: there is only a 0.4%-3.25% 
probability of these associations occurring by chance. These associations therefore require 
some explanation. Any explanation that physically links these strong X-ray sources to the 
[O III] sources does not allow the [O III] sources' PNe classifications to stand. 

Only detailed spectroscopy is likely to determine the true nature of these [O III] sources 
and the reason for their proximity to bright X-ray emitters. If (as our narrow-band flux 
estimates suggest) they are X-ray weak, young SNRs like SN 1006 or 3C 58, the existence of 
the nearby X-ray sources is still unexplained. The nearby X-ray sources cannot be ejected 
XRBs. As stated above (see §4), such a scenario requires core-collapse supernovae and 
supernova kick velocities ^ 3000 km s" 1 to explain the distances from the nebulae to the 
X-ray sources. These kick velocities would disrupt any X-ray binary system. 

Finally, perhaps an XRB/PN association could be produced within a relic of a star 
cluster. This explanation requires underlying star cluster relics in M31 massive enough not 
to dissociate completely in ~10 9 years and faint enough to be undetected in ground-based 
images. Otherwise, the best explanation for these associations appears to be the statistically 
unlikely one that all of these [O III] /X-ray pairs are spatial coincidences. 

We thank John Raymond for advising in the interpretation of the [O III] fluxes. Sup- 
port for this work was provided by NASA through grant number GO-9087 from the Space 
Telescope Science Institute and through grant number GO-3103X from the Chandra X-ray 
Center. MRG acknowledges support from NASA LTSA grant NAG5-10889. 
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Table 1. X-ray sources near cataloged PNc 



OBJ 


RA 


Dec. 


PN 


Sep (") 


Slope 


kT (kcV) a 


N H h 


X 2 /" 


Q c 


#cts 


Lx d 


rl-24 


00:42:43.211 


+41 


16:40.39 


Ford 322 


0.74 


2.1+0.2 




2.5+0.7 


17.29/19 


0.570 


227 


4.4 


rl-26 


00:42:45.095 


+41 


15:23.36 


Ford 21 


1.11 


1.79+0.03 


0.09+0.02 


2+1 


5.77/13 


0.954 


194 


6.3 


rl-2 


00:42:47.182 


+41 


16:28.56 


Ford 13 


1.24 


1.88+0.02 


0.84+0.05 


1.5+0.1 


142.5/127 


0.164 


3599 


63.8 


r3-21 


00:43:03.027 


+41 


20:41.56 


Ford 201 


1.17 


2.0+0.5 




9+4 


8.63/6 


0.195 


93 


4.0 


r3-67 


00:43:06.610 


+41 


19:13.98 


Ford 494 


2.19 


2.3+0.6 


0.05+0.01 


5+3 


1.70/3 


0.636 


73 


117 


r3-7 


00:43:21.063 


+41 


17:50.44 


Ford 209 


2.31 


2.0+0.2 




0.6 


13.44/19 


0.815 


224 


4.2 



a For sources with a significant blackbody component, the blackbody temperature is provided with a la error. 
b The absorption with la error in units of 10 21 cm" 2 . 

c The probability (based on x 1 jv) that the data represent a sample taken from a source spectrum with the model 
parameters listed. 

d The unabsorbed X-ray luminosity (0.3-7 keV) in units of 10 36 erg s -1 , assuming the spectral fit represents the 
true source spectrum. 
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Table 2. PNe Narrow Band Fluxes in units of 10 







erg 


cm 2 s~ 


-l 




PN 


Ha 


[S II] 


[O III] 


[S II]/Ha a 


[O III]/Ha b 


Ford 322 


3.4 


3.1 


3.4 


0.9 


1.3 


Ford 21 


1.1 


1.5 


11.4 


1.4 


13.4 


Ford 13 


5.9 


6.5 


9.1 


1.1 


1.9 


Ford 201 


0.9 


0.6 


3.3 


0.7 


4.7 


Ford 494 


<0.1 


<0.1 


1.0 




>10 


Ford 209 


3.4 


0.7 


11.2 


0.2 


3.6 



a A [S II] /Ha ratio greater than 0.4 is typical for shock 
heated gas. Lower ratios are typical of photo-ionized gas. 

b [0 III] /Ha ratios are corrected for absorption assum- 
ing A v = JVff/1.79 x 10 21 (Predehl & Schmitt 1995), and 
the standard extinction law. Where our fits measured 
N H > 2 x 10 21 cm" 2 , N H = 2 x 10 21 cm" 2 was assumed. 
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00 42 31.7 31.5 31.3 31.1 30.9 00 42 41.9 41.7 41.5 41.3 41.1 

RA (J2000) RA (J2000) 

Fig. 1. — Four sections of the LGS [O III] image are shown. The la error circles of the 
positions for strong X-ray sources located in 4 globular clusters are overplotted. The X-ray 
positions are clearly within the errors of the globular cluster centers, revealing the precision 
of the X-ray/optical alignment. 
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Fig. 2. — The error circles of Chandra X-ray positions of sources associated with [0 III] 
nebulae classified as PNe in M31. The la error circles are plotted on sections of an aligned, 
continuum-subtracted [0 III] image. The [0 III] nebulae show up as dark areas on the 
images, and they are labeled with their PN catalog names. 
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Fig. 3. — The results of our monte carlo tests for coincidental matches of PNs and X-ray 
sources as a function of maximum separation. Left Panel: The histogram shows the number 
of matches found in the real, aligned data. The dotted line shows the mean number of 
false associations. Right Panel: The dotted line shows the probability of 6 or more false 
associations in the randomly rotated data, and the dashed line shows the probability of 4 or 
more false associations. 



